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Abstract 

We present results from the Suzaku observations of the dwarf nova SS Cyg in 
quiescence and outburst in 2005 November. High sensitivity of the HXD PIN and 
high spectral resolution of the XIS enable us to determine plasma parameters with 
unprecedented precision. The maximum temperature of the plasma in quiescence 
20.4i2;6(stat.)±3.0(sys.) keV is significantly higher than that in outburst G.Ot^J keV. 
The elemental abundances are close to the solar ones for the medium-Z elements (Si, 
S, Ar) whereas they decline both in lighter and heavier elements, except for that of 
carbon which is 2 solar at least. The solid angle of the reflector subtending over an 
optically thin thermal plasma is Q'^/27r — 1.7 ± 0.2 (stat.) ±0.1(sys.) in quiescence. 
A 6.4 keV iron Ka line is resolved into narrow and broad components. These facts 
indicate that both the white dwarf and the accretion disk contribute to the reflection. 
We consider the standard optically thin boundary layer as the most plausible picture 
for the plasma configuration in quiescence. The solid angle of the reficctor in outburst 
Q^/2tt = O.dtli and a broad 6.4 keV iron line indicate that the reflection in outburst 
originates from the accretion disk and an equatorial accretion belt. The broad 6.4 keV 
line suggests that the optically thin thermal plasma is distributed on the accretion 



disk like solar coronae. 
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1. Introduction 

Dwarf novae (DNc) arc non-magnetic cataclysmic variables (CVs; binaries between a 
white dwarf primary and a mass-donating late-type star) which show optical outbursts typically 
with Amy = 2-5 lasting 2-20 d with intervals of ~10 d to tens of years (Warner 1995). These 
outbursts can be explained as a result of a sudden increase of mass-transfer rate within the 
accretion disk surrounding the white dwarf due to a thermal-viscous instabihty (Osaki 1974; 
Meyer & Meyer- Hofmeister 1981; Bath & Pringle 1982; Smak 1984; Cannizzo 1993; Osaki 
1996). A boundary layer (hereafter abbreviated as BL) is formed between the inner edge of 
the accretion disk and the white dwarf where matter transferred through the disk releases its 
Keplerian motion energy and settles onto the white dwarf. BL is a target of EUV and X- 
ray observations since its temperature becomes T~ 10^-10^ K. Pringle & Savonije (1979) and 
Patterson & Raymond (1985) have predicted that radiation from the BL starts to shift from 
hard X-ray to EVU when the high M front arrives at the inner edge of the disk, because BL 
becomes optically thick to its own radiation. This prediction has been verified by a number of 
multi-waveband coordinated observations (Ricketts et al. 1979; Jones & Watson 1992; Wheatley 
et al. 2003). 

The region around the inner edge of the disk is filled with a lot of intriguing but still 
unresolved issues. Within the framework of the standard accretion disk (Shakura & Syunyaev 
1973), half of the gravitational energy is released in the accretion disk, and hence, the other half 
is released in BL. The observations in extreme-ultraviolet band of VW Hyi and SS Cyg, however, 
revealed that the fractional energy radiated from BL is only < 10% of the disk luminosity 
(Mauche et al. 1991; Mauche et al. 1995). According to the classical theory, the temperature of 
BL in outburst is predicted to be 2-5 x 10^ K (Pringle & Savonije 1979), whereas the temperature 
estimated by ultraviolet and optical emission lines is constrained to a significantly lower range 
5-10x10^ K (Hoare & Drew 1991). These discrepancies may be resolved if we assume that BL 
is terminated not on the static white dwarf surface but on a rapidly rotating accretion belt 
on the equatorial surface of the white dwarf (Paczyhski 1978; Kippenhahn & Thomas 1978). 
Suggestions of the accretion belt, rotating at a speed close to the local Keplerian velocity, have 
been reported from a few DNe in outburst (Long et al. 1993; Huang et al. 1996; Sion et al. 
1996; Cheng et al. 1997; Szkody et al. 1998). Mechanism has not been understood yet to drive 
a dwarf nova oscillation (DNO) which is a highly coherent oscillation of soft X-ray and optical 
intensities with a period of 3-40 s (Robinson et al. 1978; Cordova et al. 1980; Cordova et al. 
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1984; Schoembs 1986; Marsh & Home 1998; Patterson et al. 1998; Mauche & Robinson 2001). 
Warner & Woudt (2002) try to understand DNO by assuming a magnetically driven accretion 
onto the accretion belt by enhancing a magnetic field the belt through a dynamo mechanism. 

One of the unresolved outstanding issues may be the origin of a hard X-ray optically thin 
thermal emission in outburst, since BL is believed to be optically thick. In order to identify its 
emission site, and to obtain some new insight on BL in quiescence as well, we planned to observe 
SS Cyg both in quiescence and outburst with the X-ray observatory Suzaku (Mitsuda et al. 
2007). SS Cyg is a dwarf nova in which the 1.19±O.O2M0 white dwarf and the O.7O4±O.OO2M0 
secondary star (Friend et al. 1990) are revolving in an orbit of i = 37° ± 5° (Shafter 1983) with 
a period of 6.6 hours. The distance to SS Cyg is measured to be 166 ± 12 pc using HST/FGS 
parallax (Harrison et al. 1999). SS Cyg shows an optical outburst roughly every 50 days, in 
which my changes from 12th to 8th magnitude. The optically thin to thick transition of BL has 
clearly been detected with coordinated observations of optical, EUVE, and RXTE (Wheatley 
et al. 2003). In § 2, we show an observation log and procedure of data reduction. In § 3, detail 
of our spectral analysis is explained. Owing to a high spectral resolution of the X-ray Imaging 
Spectrometer (XIS; Koyama et al. 2007) and a high sensitivity of the Hard X-ray Detector 
(HXD; Takahashi et al. 2007; Kokubun et al. 2007) over 10 keV enable us to determine spectral 
parameters of hard X-ray emission of SS Cyg with unprecedented precision. In § 4, we discuss 
on the emission site and its spatial extension both in quiescence and outburst by utilizing the 
spectral parameters, a 6.4 keV iron line parameters in particular. We summarize our results 
and discussions in § 5. 

2. Observation and data reduction 

2.1. Observations 

The Suzaku observations of SS Cygni in quiescence and outburst were carried out during 
2005 November 2 01:02-23:39(UT), and 2005 November 18 14:15(UT)-November 19 20:45(UT), 
respectively, as part of the performance verification programme (Ishida et al. 2007). In Fig. 1, 
we show an optical light curve covering our SS Cyg observations taken from the home page 
of American Association of Variable Star Observers (AAVSO)^. The outburst observation was 
performed ~two days after the optical maximum. The observation log is summarized in table 1. 

Throughout the observations, XIS was operated in the normal 5x5 and 3x3 editing modes 
during the data rate SH/H and M/L, respectively, with no window/burst options, while the 
HXD PIN was operated with the bias voltage of 500 V for all the 64 modules. We do not use 
the HXD GSO data, because we did not receive any significant signal from SS Cyg throughout. 
The source was placed at the XIS nominal position in both observations where the observation 
efficiency of all the XRTs are more than 97% of maximum throughput (Serlemitsos et al. 2007). 

^ http://www.aavso.org/ 
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Fig. 1. A V-band light curve of SS Cyg covering the Suzaku observations. The observation periods of 
quiescence and outburst are segmented with dashed vertical lines. Taken from the AAVSO home page. 



Table 1. Suzaku observation log of SS Cyg 



State 


Seq. # 


Observation Date 
[UT] 




Pointing 


Detector* 


Mode 
[ks] 


Exp.* 

[cs-i] 


Intensity^ 


Quiescence 


40006010 


2005 Nov. 02 01:02 


— 02 23:39 


XIS nom. 


FI 


Normal 


39.3 


3.308±0.005 












BI 


Normal 


39.3 


4.480±0.011 












PIN 


Normal 


27.0 


0.163±0.005 


Outburst 


40007010 


2005 Nov. 18 14:15 


— 19 20:45 


XIS nom. 


FI 


Normal 


56.0 


1.580±0.003 












BI 


Normal 


56.0 


2.427±0.007 












PIN 


Normal 


47.9 


0.028±0.003 



# FI: Frontsidc-IUuminated CCD (XISO, 2, 3), BI: Backside-Illuminated CCD (XISl), PIN; HXD PIN detector. 



* Exposure time after data screening. 

t Intensity in the 0.4-10 keV (FI and BI CCDs) and 12-40 keV bands (PIN) with a Icr statistical errors. The systematic error of the PIN 
background (5% of the NXB) is 0.024 c s-Mn this band. 

As already known from previous observations (e.g.Wheatley et al. 2003) the hard X-ray flux 
was smaller in outburst than in quiescence. 

2.2. Data screening 

In data reduction, we use event files that are created with the pipeline processing software 
(revision 1.2.2.3), and the analysis software package HEASOFT (version 6.1.3). For the XIS, 
we selected the events with grades 0, 2, 3, 4, and 6. We do not use bad pixels or columns where 
charge transfer efficiency is too low to detect/transfer X-ray signals. We also discarded the 
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data taken during time intervals of maneuvers, low data rate, passage of SAA (South Atlantic 
Anomaly), the field of view being occulted by the earth or watching bright earth rim, and the 
pointing position being away from the source by more than 1^5. Finally, the cleaned event files 
are created by removing hot /flickering pixels. In addition to the criteria described above, we 
accept the data without telemetry saturation. After these data selection, we combined the data 
of 3x3 and 5x5 editing modes before starting analysis. Resultant exposure time is 39.3 ks and 
56.0 ks in quiescence and outburst, respectively. For creating a cleaned event file for the HXD 
PIN, on the other hand, we further removed the time interval during cut-off rigidity is less than 
6 GeV c-^ 

In extracting photons from the source with the XIS, we take a circular integration region 
with a radius of 250 pixels (4^34) centered at the source, which includes more than 96 % of 
the X-ray events, while an annulus with an outer radius of 432 pixels (7'50) around the source 
region is adopted for the background-integration region. In view of statistics, we took the area 
of the background-integration region twice as large as that of the source-integration region. For 
the HXD-PIN, on the other hand, we accumulate a non- X-ray background (NXB) spectrum 
from a simulated NXB event file^ published by the HXD team. In addition to the NXB, we 
need to consider the cosmic X-ray background (CXB), for which wc adopt the empirical model 
spectrum constructed on the basis of the HE AO observations (Boldt 1987), 

/cxb(B) = 9.0x(— ) xexp(-— ) erg cm- str"' keV"', (1) 

and create a CXB spectrum from this model by the f akeit command in the spectral fit software 
XSPEC (Arnaud 1996) with an exposure time of 1 Ms. In doing this, we adopt the PIN flat 
sky response. The background spectrum for the PIN is created by combining the NXB and 
CXB spectra using the mathpha command in the FTOOLS package. Since the CXB level is 
~5% of the NXB level, we ignore sky-to-sky variation of the CXB. The counting rates listed in 
table 1 are those after all the data screening described in this section are applied. 

2.3. Light Curves 

We show energy-resolved and background-subtracted light curves of SS Cyg in quies- 
cence and outburst in Fig. 2. The source and background photons are accumulated using the 
integration regions explained in § 2.2. The data from the four XIS modules are combined to 
generate the light curves in the three energy bands 0.2-0.4 keV, 0.4-2.0 keV, and 2.0-10.0 keV. 
The light curve of the HXD (PIN) is created in the energy range 12-40 keV where the sensi- 
tivity of HXD becomes the maximum. The mean counting rates are summarized in table 1. In 
outburst, the average intensity in the 12-40 keV band is 5.9% of the NXB intensity of the PIN, 
whereas there remains 5% systematic error in the NXB model intensity. Hence, the detection 
of the PIN in outburst is marginal. In the energy bands above 0.4 keV, the source is brighter 

^ ftp://ftp.darts.jaxa.jp/pub/suzaku/verl. 2/ 
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Fig. 2. Background-subtracted light curves of SS Cyg in quiescence (left) and outburst (right) in the 
bands (a) 12-40 keV, (b) 2 10 keV, (c) 0.4-2 keV, and (d) 0.2-0.4 keV. Data from all the XIS modules 
are combined to make the light curves in the panels (b) through (d). The bin size of the XIS is 256 sec 
both in quiescence and outburst whereas that of the HXD PIN is 1024 sec in quiescence and 4096 sec in 
outburst. The origin of time is MJD 57676 01:23:49 and MJD 57692 14:56:10, respectively. 

in quiescence than in outburst, as demonstrated by Wheatley et al. (2003), and is more vari- 
able. The flux below 0.4 keV is, on the other hand, higher in outburst. Moreover, the source 
obviously declines throughout the observation, whereas the intensity is nearly constant in the 
higher energy bands. These facts suggest that Suzaku detected high energy end of the emission 
from the optically thick BL which appears in outburst (Tylenda 1977; Pringle 1977; Pringle & 
Savonije 1979; Patterson & Raymond 1985; Wheatley et al. 2003). 

2.4- Average spectra 

In Fig. 3, we show averaged spectra of SS Cyg from the XIS and the HXD-PIN in quies- 
cence and outburst. The source and background photons are accumulated using the integration 
regions explained in § 2.2. The data from XISO, 2, and 3 are combined into a single FI-CCD 
spectrum. The BI-CCD spectrum originates solely from the XISl data. 

SS Cyg is detected at least up to ~30 keV with the HXD-PIN in quiescence, whereas 
the detection of the PIN in outburst is marginal as noted in § 2.3. From the inset, iron Ka 
emission lines are clearly resolved into 6.4, 6.7, and 7.0 keV components. The 6.4 keV line 
indicates reflection of the hard X-ray emission from the white dwarf and possibly from the 
accretion disk, as pointed out by Done & Osborne (1997). Note that the 6.4 keV line is broad 
in outburst. Except for the iron emission lines, only weak signs of H-like Ka lines are visible 
from the other elements in quiescence. The outburst spectra, on the other hand, are softer 
than those in quiescence, and are characterized by H-like and He-like Ka emission lines from 
nitrogen to iron (Okada et al. 2008). From the inset, the He-like iron Ka line at 6.7 keV is 
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Fig. 3. Averaged spectra of SS Cyg in quiescence and outburst. The source and background integration 
regions are explained in § 2.2. Tlie data from XISO, 2, and 3 are combined into a single FI-CCD spectrum 
whereas BI-CCD spectrum originates from XISl. The insets are the blowups in the 5-8 keV band showing 
profiles of iron Ka emission lines. 

much stronger than the H-hke hne at 7.0 keV, in contrast to the quiescence spectra. These 
facts indicate that the plasma has a temperature distribution, and that the average plasma 
temperature is significantly lower in outburst than in quiescence. The decaying soft component 
in the light curves (Fig. 2) appears as an excess soft emission below ~0.3 keV in the BI-CCD 
spectrum in outburst. 

3. Analysis and Results 

3.1. Spectral components and their models 

In Fig. 4, we show schematic view of spectral ingredients of SS Cyg in the 0.2-40 keV 
band. This picture is originally proposed by Done & Osborne (1997) on the basis of their 
analysis of the Ginga and ASCA data of SS Cyg. The hard spectrum of SS Cyg is primarily 
composed of an optically thin thermal plasma emission with a temperature distribution (§ 2.4), 
for which Done & Osborne (1997) introduced a power-law type differential emission measure 
(DEM) model as 

d(EM) « d{\ogT) oc ' dT, (2) 

V max / V max / 

where Tmax is the maximum temperature of the plasma. The model they used, named cevmkl 
in XSPEC, is an optically thin thermal plasma emission model (mekal in XSPEC, Mewe et 
al. 1985; Mewe et al. 1986; Liedahl et al. 1995; Kaastra et al. 1996) convoluted with this 
temperature distribution. Since it can successfully represent ~30 spectra of the dwarf novae 
(Baskill et al. 2005) observed with ASCA, we adopt this model as well in this paper. In addition 
to this main component, its reflection from the white dwarf (and the accretion disk) occupies 
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Fig. 4. Schematic view of the spectral component of SS Cyg in the 0.2-40 keV band. 

significant fraction of the observed X-ray fiux, as required from detection of a fiuorescent iron 
Ka line at 6.4 keV. To represent the refiection, we adopt the model reflect (Magdziarz & 
Zdziarski 1995), which is a convolution-type model describing refiectivity of neutral material. 
In summary, the X-ray spectra of SS Cyg are composed of the multi-temperature optically thin 
thermal plasma component, its refiection from the white dwarf (and the accretion disk), and 
the fiuorescent iron emission line at 6.4 keV. 

3.2. Evaluation of plasma parameters 

Our purpose is to constrain the geometry of the optically thin thermal plasma in SS Cyg. 
In doing this, wc can utilize the equivalent width (EW) and the profile of the 6.4 keV emission 
line. The EW reflects a covering fraction of a reflector viewed from the plasma (Makishima 
1986). The proflle contains information of motion of the reflector. From continuum analysis, 
the covering fraction of the reflector can independently be evaluated. 

To utilize these methods, however, it is essential to know a priori the iron abundance 
Zpe. For its estimation, we generally make use of the relative intensities and the EWs of the 
He-like and H-like iron Ka emission lines at 6.7 keV and 7.0 keV, respectively. In deriving 
from these quantities, we need to know the plasma emission parameters such as Tmax and 
a through spectral fitting of the cevmkl model to the observed data. One may suspect that 
it is enough to fit the spectrum containing the He-like and H-like lines locally with a single 
temperature optically thin thermal plasma model. This does not work properly, however, since 
the temperature of the plasma in SS Cyg is distributed in such a wide range, from Tmax — 20 keV 
(see below) down to ~0.1 keV indicated by the nitrogen and oxygen lines, that plasma emission 
components whose temperatures are too high enhance only the continuum level and dilute the 
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iron emission lines, thereby resulting in a lower abundance estimation (Done & Osborne 1997). 

To know the plasma emission parameters Tmax and a, however, we need to know the 
metal abundances, the iron abundance (^Pe) in particular and the oxygen abundance (Zq) as 
well, and the covering fraction of the reflector (Q), because they couple with each other in 
the spectrum above 7 keV. At this stage, we have found that the plasma emission parameters 
(^max and a) and the reflection parameters (Q, Zpe and Zq) depend on each other, and need to 
determine them in a self-consistent manner. Accordingly, we carried out a combined spectral 
fit among the quiescence and outburst spectra in selected bands crucial to determine T^ax, 
Q, Z-pc, and Zq- They are the quiescence spectra in the 4.2-40 keV band (sensitive to T^^^, 
a^, Zpc, and fl^), the outburst spectra in the 4.2-10 keV (sensitive to T^^^, a^, Zpc, and O*-*) 
and that in 0.55-0.72 keV bands (sensitive to Zq and q;°), where the superfixes 'Q' and 'O' 
indicate quiescence and outburst, respectively. We ignore the intermediate 0.72-4.2 keV band 
in outburst, because (1) this band is of no use to evaluate the parameters listed above. Including 
this band may rather introduce unexpected systematic errors in the parameters of interest, (2) 
this band includes the emission lines from Nc to Ca whose abundances arc uncertain until we 
determine and a'^, and (3) as will be found later in table 2 and 3, values of a*-* are different 
for different elements. The entire outburst spectrum can not be fit with a single a cevmkl model. 
We set Tmax, ct, and Q free to vary in quiescence and outburst separately, because the geometry 
of the plasma, and hence its temperature distribution may be different, whereas we set Z^e and 
Zq common between quiescence and outburst. Note that we do not use the PIN spectrum in 
outburst, because the detection is marginal due to the NXB uncertainty (§ 2.3 and 2.4). The 
result of the fits is shown in Fig. 5, and its best-fit parameters are summarized in table 2. The 



Fig. 5. Result of the simultaneous fits of the cevmkl plus reflection model to the 4.2-40 keV band in 

quiescence, the 4.2-10 keV and 0.55 0.72 keV bands in outburst. Black, red, and green colors are used 
for the data and the model of the FI-CCD, BI-CCD, and PIN, respectively, except for the outburst 
0.55-0.72 keV spectrum where only the BI-CCD data are utilized. 

fits are marginally acceptable at the 90% confidence level. The parameters of the plasma and 
the refiector are determined consistently between quiescence and outburst with unprecedented 
precision. We adopt Anders & Grevesse (1989) as the solar abundances of the metals. The DEM 
power a is varied independently between the iron and oxygen bands in the outburst spectra, 
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Table 2. Best-fit parameters of the simultaneous fit of the cevmkl plus reflection model to the quiescence and outburst 
spectra. 



Phase Quiescence Outburst 

Energy band [keV] 4.2-40 4.2-10 0.55-0.72 

T„.ax* [keV] 20.4t^-0(±3.o) 6.0t?;^ 

at 0.7t°:?(±0.1) 5.8t?-i -O.SetH; 

f2/27rt 1.7 ±0.2 (±0.1) O.Qtoi 

§ 0.37tn3(±0-01) 

^o" 0.46tn3(±0-01) 

(d.o.f.) 728.2 (657) 



JNote — Errors in parentheses are systematic errors associated 

with the NXB uncertainty of the HXD-PIN. 

* Maximum temperature of the optically thin thermal plasma. 

t Power of DEM as d(EM) oc (T/Tmax)«-i(iT. 

t Solid angle (covering fraction) of the reflector viewed from the 

plasma. 

§ Iron abundance in a unit of solar. Constrained to be common 
among the three energy bands. The solar abundance table of 
Anders & Grevesse (1989) is adopted, in which [Fe/H] = 4.68 x 

W-''. 

I Oxygen abundance in a unit of solar. Constrained to be common 
among the three energy bands. The solar abundance table of 
Anders & Grevesse (1989) is adopted, in which [O/H] = 8.51 x 
10-*. 

because the cooling functions are significantly different between these two bands (Gehrels & 
Williams 1993). In fact, the resultant a values are significantly different between O (~ —0.4) 
and Fe (~6). Considering this large difference, we have checked mutual consistency of the model 
spectra in the two energy bands in outburst. The result is shown in Fig. 6. The extrapolation of 
the best-fit model in one band does not exceed the observed flux in the other band. Hence, the 
flts to the two energy bands in outburst are revealed to be mutually consistent. The maximum 
temperature of the plasma in quiescence (T^^^x — '^^■^-ts keV) is significantly higher than that 
in outburst (T^^^ = Q-O-ii keV). Our T^^ is consistent with that of Done & Osborne (1997) 
(= 2115^7 keV). Their T^g^ (= Q-Qtfj keV) seems slightly higher than ours, though this may 
be due to different outburst phases at which the observations were carried out. The optically 
thin thermal plasma is covered with the refiector more in quiescence (Q^/27r = 1.7 ± 0.2) than 
in outburst {OP /2t^ — 0.9lo;4)- The covering fractions are not constrained by Done & Osborne 
(1997) very weU (QQ/27r = O.G^oie and VP /2t: = 2.2t\l). The iron and oxygen abundances are 
both sub solar (Zpe = 0.371^3^© and Zq = 0.461^3^0)- 

Finally, we have considered the effects of the 5% systematic PIN NXB error. We have 
made background spectra of the PIN with an intensity being enhanced or reduced by 5%, 
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Fig. 6. The outburst spectra and the models used to determine the parameters of the plasma emission 
and the reflection. Only the data of the XISl are shown for clarity in red. The green and the blue curves 
show the best-fit models in the 0.55-0.72 keV and 4.2-40 keV bands, respectively, whose parameters are 
summarized in table 5. 

subtracted them from the quiescence PIN spectrum, and repeated the combined spectral fit 
described above. The resultant systematic errors are summarized in parentheses in table 2. 
We have found T^^^x accompanied by the largest fractional systematic error of ±3.0 keV. 
Fortunately, those of the other parameters are smaller than the statistical ones. Since we do 
not use the outburst PIN data, there appears no systematic error in T^^^, a^, and VP /2t{. 

3.3. Abundance of the other elements 

We can identify He-like and H-like Ka emission lines from nitrogen to iron in the outburst 
spectra shown in Fig. 3. Utilizing the intensities of these lines, we can estimate the abundances 
of corresponding elements. The line intensities depend upon Tj^g^^, and the abundance. 
Among them. We have constrained T^^^ well by the combined fit analysis described in § 3.2. 
With r°g^ being fixed at 6.0 keV (table 2), we have evaluated the abundances by fitting the 
cevmkl model to the Ka lines of each element separately. In principle, (see eq. (2)) is 
determined by the relative intensities of the He-like and H-like Ka lines, and the abundance 
is constrained by and their EWs. The results of the fits are shown in Fig. 7 and the best- 
fit abundances as well as are summarized in table 3. The spectral fit is started from 
a higher energy band. This is because the line emissions from lighter elements do not affect 
the determination of the abundances of heavier elements, whereas L-shell emission from the 
heavier elements is possible to affect the fit to the lighter element lines. As the abundances 
of the heavier elements have been measured, those of the lighter elements are determined one 
after the other by fixing the abundances of the heavier elements at their best-fit values. The 
fit to Ca and Ar is carried out simultaneously, and so is to N and C. Note that only BI-CCD 
is used for the fit to N and C. 

In order to estimate the abundances of N and C, we have to take into account several 
possible systematic effects such as the low energy response of the XIS, the uncertainty of the 
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Fig. 7. Local cevmkl fits to the He-like and H-like Ka emission lines of Ca, Ar, S, Si, Mg, Ne, N, and C. 
The upper panel shows the data (crosses) and the models (histograms), while the lower panel shows the 
residuals in a unit of the standard deviation. In order to show detection of the C emission lines, we have 
also plotted a cevmkl model with the C abundance being set equal to zero as a dotted histogram in panel 

(f)- 
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Table 3. Summary of the abundances of iron, calcium, argon, sulfur, silicon, magnesium, neon, oxygen, and carbon. 



Element 


Energy Band 
(keV) 


(10i9cm-2) 


ro t 

max 

(keV) 




Abundance 
(Zo) 


(d.o.f.) 




Fe^ 


4.2-10 


5.0f 


6.ot;j 




o.37t°:[!J 


728.3 (657) 


1.11 


Ca 


2.8-5.0 


5.0^ 


6. of 


9 9+0.5 


"•'•^-0.36 


263.7 (265) 


1.00 


Ar 










q3+"--^4 

'-'■^"->-0.32 






S 


2.15-3.0 


5.0^ 


6.0'' 




0.88±0.09 


242.1 (171) 


1.42 


Si 


1.6-2.3 


5.0^ 


6.0*' 




0.90 ±0.07 


389.0 (244) 


1.59 


Mg 


1.2-1.7 


5.0^ 


6.0^ 


0.44 ±0.06 


0.80 ±0.06 


316.2 (267) 


1.19 


Ne 


0.8-1.2 


5.0^ 


6.0^ 


0.31 ±0.02 


0.49 ±0.04 


429.8 (217) 


1.99 


fit with the BI spectrum only 


0^ 


0.55-0.72 


5.0^ 






46+"-°'' 


728.3 (657) 


1.11 


N 


0.3-0.52 


3.5f 


6.0' 




U-*U_Q 4Q 


44.7 (49) 


0.91 


C 










47+1.60 

-0.99 






N 


0.3-0.52 


5.0f 


e.of 


64+° " 




44.6 (49) 


0.91 


C 










•^•^^-0.90 






N 


0.3-0.52 


7.9^ 


6. of 


56+°-^° 


"•■■-"^-o.is 


44.6 (49) 


0.91 


C 










9 7C+1.14 







Note. — 'f implies the parameter is lixed. 

t Fixed at 6.0 keV obtained from the simultaneous fit (table 2 in § 3.2). 

t Fixed at 5.0 x 10^^ cm~^ obtained from the Chandra LETG observation (Mauche 2004). 

' The same values as summarized in table 2. 



contamination accumulating on the optical blocking filters above the CCD chips, and systematic 
error of the hydrogen column density to SS Cyg. Among them, the low energy response and 
the contamination are investigated by ourselves with the PKS 2155-304 data taken during Nov. 
30 - Dec. 2, 2005, which date is close enough to our SS Cyg observations. As summarized 
in Appendix 1, we need to apply additional carbon K-edge with r = 0.88 ± 0.05 and extra 
iVn = 8.2 ±0.7 X 10^^ cm~^. By applying these corrections to the model, we can utilize the BI- 
CCD data down to 0.23 keV. In order to estimate the A^h systematic error, on the other hand, 
we have checked variation of the best-fit parameters with an A^h of 3.5, 5.0, and 7.9 x 10^^ cm~^. 
These values are selected based on A^h = 5.0^^5 x 10^^ cm~^ obtained from the Chandra LETG 
observation of SS Cyg in outburst (Mauche 2004). Since the H-likc Ka line of C and the He-like 
Ka line of N cannot be resolved with the BI-CCD energy resolution, we employ the BI-CCD 
spectrum in the 0.3-0.52 keV band covering the emission lines of both C and N, and fit cevmkl 
model. 

We remark that not all the fits shown in Fig. 7 are acceptable. The fit to Si is especially 
poor. As evident from the fit residuals inconsistent between the BI and FI-CCDs at around 
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the Si K-edge (1.7-1.9 keV), this is mainly due to systematic error associated with energy 
response cahbration. Nevertheless, the best-fit parameters summarized in table 3 indicate that 
the abundances of the medium-Z elements (Si, S, and Ar) are close to solar ones, whereas they 
are reduced for both lighter and heavier elements. The exception is the carbon abundance 
whose lower limit is 2.0Zq even if we consider the systematic error of A^h- Since the energy 
resolution of the XIS in the C line energy band is not so good as in higher energy band, we 
have drawn the cevmkl model with the C abundance being set equal to zero in Fig. 7(f) as the 
dotted histogram. The data excess associated with the C emission is clear. 

Given the spectral parameters, we can calculate X-ray fluxes of SS Cyg, which is 19.5 x 
10~^^ erg cm~^ s~^ and 8.1 x 10~^^ erg cm~^ s~^ in the 0.4-10 keV band in quiescence and 
outburst, respectively. Assuming the distance of 166 pc (Harrison et al. 1999), we obtain the 
luminosities of 6.4 x 10^^ erg s~^ and 2.7 x 10^^ erg s~^, respectively, in this energy band. 

3.4- Iron emission lines 

3.4-1- Quiescence 

In order to constrain the geometry of the optically thin thermal plasma in quiescence, 
we have attempted to evaluate the parameters of the 6.4 keV line in detail. We fit a model 
composed of a power-law continuum and three narrow Gaussian lines at 6.4 keV, 6.7 keV and 
7.0 keV to the XIS spectra in the 5-9 keV band. The result of the fit is shown in Fig. 8 and the 
best-fit parameters are summarized in table 4. There remains a slight excess emission in the 



Fig. 8. Fits of a power law plus 3 Gaussians (left) and 4 Gaussians (right) to the 5-9 keV XIS spectra in 
quiescence. 

low energy tail of the 6.4 kcV line. We thus have added a broad Gaussian to the 6.4 keV line. 
The results are also summarized in table 4 and Fig. 8. The fit is improved with a (d.o.f) value 
from 261.7 (233) to 251.5 (231). The chance probability for this to occur is 0.009, which implies 
adding a broad Gaussian is significant at ~99% confidence level. The resultant equivalent width 
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Table 4. Best-fit parameters of the iron Ka emission lines with the 5-9 keV XIS spectra in quiescence. 



Component 


Parameter 


3 lines 


4 lines 


Continuum 




l-60_o.o2 


1 GO+0.02 
i.OZ_o 03 




AT t 
JVC' 


y.io ± U.U/ 


y.ou ± U.Us 


H-like line 


El,* (keV) 


A/-'o+0 003 

6-968+0.008 


A/-'/i+0 006 

6-966^0.008 




AT t 


c .)n+0.39 
0.>3U_o.41 


c •)Q+0-38 
O.oy_o.42 




a [KB V j 


U (nxeaj 


U (nxeaj 




bjVv^ (eVj 


1 on+9 


-'■'^"^-10 


He-like line 


El,* (keV) 


6.682^^0^ 


6.684inS^ 




Ni,i 


7 O1+0.41 
'•°-'--0.42 


7 fi7+0.36 
' -0.48 




a (keV) 


(fixed) 


(fixed) 




EW§ (eV) 


181 ±10 




Fluorescence line parameters 


Narrow line 


El,* (keV) 


fi 407+0-008 


fi 404+0-012 




Ni,i 


*-U^-0.36 


2 4fi+0-9i 




a (keV) 


(fixed) 


(fixed) 




EW§ (eV) 


86±8 


53lf 


Broad line 


Ni,i 




9 qo+0.98 
Z.OO_o 50 




a (keV) 




n 1 1 9+0.055 
U.llz_o 031 




EW§ (eV) 








(d.o.f) 


261.7 (233) 


251.5 (231) 




xl 


1.12 


1.09 



' Line central energy, i'he -value is constrained to be the same 

between the narrow and broad 6.4 keV components. 

t Continuum normalization in a unit of 10~^ cm~^ keV~^ s"'^ at 



I keV. 

t Line normalization (intensity) in a unit of 10~^ pho- 
tons cm~^ s~^. 
§ Equivalent width. 

II Photon index. 

of the narrow and broad lines is EW (narrow) = SSlg^ eV and EW (broad) = Slln eV. 
34.2. Outburst 

In outburst, tiie 6.4 keV line is obviously broad (Fig. 3). We therefore adopt a broad 
Gaussian for modehng the 6.4 keV emission line. As in the case of the fit to the quiescence 
spectra (§ 3.4.1), we have tried to fit a model composed of a power law and three Gaussians to 
the XIS spectra in the 5-9 keV band. The result is shown in Fig. 9, and the best-fit parameters 
are summarized in table 5. The equivalent width of the 6.4 keV line is as large as ~210 eV. 
The energy width of the fine a — 0.12 keV is nearly the same as that of the broad component 
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Table 5. Best-fit pajrameters of the iron Ka emission lines with the 5-9 keV XIS spectra in outburst. 



Component Parameter 



Continuum 


rii 




o if;+0.09 
^•^"^-0.05 








R 77+0-65 
D. / /_o.l2 


H-like Ka 


El* 


(keV) 


o.yoo_Q 




N^i 




0.82±0.17 




a (ki 


eV) 


(fixed) 




EW^ 


' (eV) 


80 ±17 


He-like Ka 


-El* 


(keV) 


6.6821°;;]°^ 








4.3710:^9 




a (b 


eV) 


(fixed) 






i (eV) 


390^2? 


Neutral Ka 


El* 


(keV) 


fi 448+0 024 

0.4^0_Q Q25 




iVL* 




2 Qn+O-95 
"^•^^-0.38 




£7 (b 


eV) 


'-'•-'^^-0.03 




EW^ 


i (eV) 


208t^f 



(d.o.f) 198.1 (155) 

Xl 1.28 
* Line central energy. 'I'he value is constrained to be the same 

between the narrow and broad 6.4 keV components. 

t Continuum normalization in a unit of 10~^ cm~^ keV~^ at 

I kcV. 

f Line normalization (intensity) in a unit of 10~^ pho- 
tons cm~^ s~^. 
§ Equivalent width. 

II Photon index. 

in quiescence. 

In quiescence, the narrow component dominates the 6.4 keV hue. In order to see if there 
is a narrow component also in the outburst spectra, we have implemented another narrow 
Gaussian (a = 0) into the model at the same central energy of the broad 6.4 keV component. 
The fit with this model, however, only provides an upper limit to the narrow 6.4 keV line with 
an equivalent width of < 39 eV, or < 20% of the broad component. 

4. Discussion 

4..!. Location of the optically thin thermal plasma in quiescence 

4.1.1. Geometry of the plasma 

We have shown in § 3.2 that the optically thin thermal plasma in quiescence is covered 
with the reflector with a solid angle of n^/27r — 1.7 ±0.2. Such a high value of il'^/27r is 
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Fig. 9. Fits of a power law plus 3 Gaussians to the 5-9 keV XIS spectra in outburst. 

achieved only in a limited number of high-mass X-ray binaries (e.g. Watanabe et al. 2003) that 
are surrounded by matter as thick as ~10^'^ H cm^^. The thickness of matter to SS Cyg, on 
the other hand, only amounts to A^h ^ 10^*^ cm^^ or r ^ 10"^ for Thomson scattering, which 
is too transparent to reflect X-rays efficiently. Candidates of the reflection site are therefore 
limited to the white dwarf and the accretion disk. The observed f2^/27r, however, requires a 
special geometry, because even an infinite slab can only subtend a solid angle of r2/27r = 1 over 
an X-ray source above it. The plasma should be compact enough compared to the size of the 
white dwarf, and should be located very close to the white dwarf and the accretion disk. We 
consider the standard optically thin BL (Patterson & Raymond 1985) as the most plausible 
configuration of the optically thin thermal plasma in quiescence (see also Fig. 10). In § 3.4.1, 



Fig. 10. Geometry of the BL in quiescence deduced from the Suzaku observation. 
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we have shown that the 6.4 keV iron hne profile favors the broad component in addition to 
the narrow component, which can naturally be interpreted as originating from the accretion 
disk and the white dwarf, respectively, via fluorescence. The width of the broad component 
a — O.llloios keV corresponds to a line-of-sight velocity dispersion of 5300li5oo km s~^, which 
is consistent with the line-of-sight Keplerian velocity amplitude of the accretion disk just on 
the 1.19M0 white dwarf (vRsini = 3800±400 km s-\ where we set i = 37° ± 5°). 
4.I.2. Size of the boundary layer 

Now that we know the abundance of iron O.37lo!o3^0 (§ 3.2), we can estimate the height 
of the boundary layer in quiescence on the basis of the standard BL geometry displayed in 
Fig. 10. George & Fabian (1991) theoretically calculated the equivalent width of the fluorescent 
iron Ka line due to a point source being located above an infinite slab (i.e, VL = 2ti). The 
equivalent width depends on the inclination angle i between the observer's line of sight and 
the normal of the slab, a photon index of the spectrum of the illuminating point source, and 
an iron abundance of the slab. In applying their calculation to the observed narrow 6.4 keV 
component, we have adopted the following parameters; 

1. We set the inclination angle i = 60°, which is the average over the visible hemisphere of 
the white dwarf surface (Done & Osborne 1997). 

2. The observed continuum spectra above 5 keV are represented with a power law with a 
photon index of ~1.6 (table 4). This photon index is, however, affected by the reflected 
continuum. In order to correct this, we have retried a flt with a power law plus its reflection 
with Vfi 1271 and Z^c flxed at the values listed in table 2. The resultant photon index is 
r = 1.8, which we adopt for the power law incident on the white dwarf surface. 

With this parameter set, we have obtained an expected equivalent width of the narrow 6.4 keV 
iron line to be 110 cV, according to Fig. 14 of George & Fabian (1991). Note here that this 
figure is drawn under the condition of the solar abundance with a composition of [Fe/H] = 
3.2 X lO^'""' whereas our case is 0.37^0 under the condition of [Fe/H] = 4.68 x 10~^ (Anders & 
Grcvcssc 1989). After correcting this abundance difference using their Fig. 16, the expected 
equivalent width of the narrow 6.4 keV line is 



where Qwd is the solid angle of the white dwarf viewed from the BL plasma. Equating this to 
the observed EW = ^?>tf (table 4) , we obtain Qwd /27r = 0.66lo:?i ■ By assuming that the plasma 
is point-like and is located at a height h above the white dwarf, we can directly link Qwd/Stt 
with h. From a simple geometrical consideration, we obtain /i = 0.060l;o;o56 or h < 0.12i?wD- 
Note that this number is similar to the thickness of the BL estimated in the echpsing dwarf 
nova HT Gas h < 0.15i?wD (Mukai et al. 1997). The total sohd angle Q^/27r can be evaluated 
also from the 6.4 keV line by comparing the total equivalent width of EW = 104 to eq. (3) as 
Jl/27r = 1.3lo;2- This is consistent with Q*^/27r = 1.7 ±0.2 estimated by the continuum spectra. 




(3) 
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We therefore conclude that the Suzaku results in quiescence favor the geometry of the standard 
optically thin BL (Patterson & Raymond 1985) and its scale height from the white dwarf is 
h < 0.12i?wD- 

4.2. Location of the thin thermal plasma in outburst 

4-2.1. Geometry of the plasma 

As presented in § 3.4.2, the 6.4 keV iron line in outburst is broad with a Gaussian a 
of 0.12 keV (table 5). Since this width can be interpreted as the line-of-sight Kepler velocity 
amplitude of the accretion disk on the white dwarf (§ 4.1.1), the main reflector in outburst is 
likely to be the accretion disk. From the upper limit of the narrow 6.4 keV line, contribution 
from the static white dwarf surface evident in quiescence is negligible with an upper limit of 
20% of the disk contribution. Moreover, the solid angle of the reflector Q°/27r = 0.9^0 4 is 
consistent with an infinite plane, reminiscent of the disk, although the error is large. All these 
facts strongly suggest that the plasma is located above the disk like coronae with their height 
small enough compared with the disk radius, as shown in Fig. 11. This picture is supported 



Fig. 11. Geometry of the optically thin thermal plasma in outburst. The energy width of the 6.4 keV 
line, the small contribution from the white dwarf surface to the EW, and 0^/277 = 0.9^q4 strongly suggest 
that the plasma is located on the disk like coronae. 

also by the Chandra HETG observations (Okada et al. 2008), in which the plasma emission 
hues are very broad with their Gaussian a as large as ~2000 km s~^. This is much broader 
than the thermal velocity, and can be attributed to the plasma being anchored to the disk by, 
for example, magnetic field and corotating with the disk. A similar disk corona geometry is 
suggested to explain the X-ray spectrum and the energy width of the emission lines in WZ Sge 
in outburst (Wheatley & Mauche 2005). 

Note that the absence of the narrow 6.4 keV component can alternatively be understood 
by invoking the accretion belt (Paczyhski 1978; Kippenhahn & Thomas 1978) which is an 
equatorial part of the white dwarf atmosphere being accelerated by the accretion torque and 
rotating nearly at the local Keplerian velocity (fxsini ~4000 km s~^ in SS Cyg, § 4.1.1). It is 
suggested from observations that the belt covers significant fraction of the white dwarf surface 
in outburst (Long et al. 1993; Cheng et al. 1997; Szkody et al. 1998). If so, the 6.4 keV hue 
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from the white dwarf should also be broad. Our data of the 6.4 keV line in outburst can be 

compatible with the existence of the accretion belt. 

4-2.2. A possible solution to the large EW of the 6.4 keV line 

As explained in § 4.2.1, the width of the 6.4 keV hne can be attributed to the Keplerian 
motion of the disk (and the accretion belt). The observed equivalent width £'1^ = 208^27 
(table 5), however, can not be explained solely by simple fluorescence due to continuum illu- 
mination. Assuming the inclination angle of the disk i — 37° (Shafter 1983), the photon index 
of the incident spectrum F = 2.3 (correction is made to the value 2.16 in table 5, see § 4.1.2), 
the iron abundance of 0.37^©, we expect EW ~ 110 eV at most according to George & Fabian 
(1991), even if we assume the upper limit of the allowed range of the solid angle Q'~'/2tt = 1.4, 
which may be achieved by taking the accretion belt into account. Since the thickness of matter 
surrounding SS Cyg is thin, only of order A^h ^ 10^° cm~^ (§ 4.1.1), there is no other source 
that can produce 6.4 keV line photons via fluorescence. The large observed equivalent width 
therefore requires some other mechanism to enhance the intensity around 6.4 keV. 

One possible solution is to invoke a Compton shoulder which appears in a lower energy 
side of a main line due to Compton down scattering. Such a shoulder is resolved with the 
Chandra HETG in a wide variety of sources such as a HMXB (Watanabe et al. 2003), a ULX 
(Bianchi et al. 2002), and an AGN (Kaspi et al. 2002). Note that the seed photon of the shoulder 
in these objects is the 6.4 keV fluorescence hne born in Compton-thick clouds, whereas in the 
current case, we consider the He-like Ka line at 6.67 keV, because it is very strong with an 
equivalent width of 390 eV (table 5). The Compton shoulder of the 6.67 keV line extends 
down to Ao -I- 2Ac in wavelength, corresponding to 6.50 keV in energy, in the case of the back 
scattering, where Aq and Ac are the wavelength of the 6.67 keV photon and the Compton 
wavelength (= h/niec), respectively. The low energy end of the shoulder (6.50 keV) is close 
enough to be merged into the broad 6.4 keV line with a CCD resolution. 

In order to evaluate contribution of the 6.67 keV Compton shoulder to the observed EW 
of the 6.4 keV line, we consider a simple geometrical model, in which a point source of 6.67 keV 
photons is centered at a height h above a circular slab with a radius of it!out- We adopt an 
inclination of the slab to be i = 37° (Shafter 1983). Then, the energy of a scattered photon 
arriving at the observer from any part of the slab can be calculated, because the energy of 
a scattered photon is uniquely determined by the scattering angle. In calculating the photon 
energy, we can ignore the Kepler motion of the slab because the situation is in a non-relativistic 
domain. Assuming that the intrinsic line emission is isotropic, we have calculated the spectrum 
of the Compton-scattered photons by integrating the emissions from the entire slab. Note 
that the intensity of the spectrum is affected by the elemental abundance of the slab through 
probability of the scattering (vs. photoelectric absorption), and absorption of the scattered 
photons along the path to escape from the slab. We adopt the abundances which we have 
already measured (§ 3.3). As a result, the cross section of the photoelectric absorption at 
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6.67 keV is obtained to be 1.15 x 10~^^ cm~^, which is roughly twice as large as that of the 
Compton scattering. The resultant iron 6.67 keV line spectrum in the case oi h — 0.1it!out is 
shown in Fig. 12. The low energy end of the shoulder appears at 6.50 keV as expected, and 



Fig. 12. Example of the line profile after a single Compton scattering by electrons in the geometrically 

thin non-rclativistic disk. The energy of the source photons is 6.67 kcV, which wc consider a He-like iron 
Ka line. The inclination angle of the disk is 37°. The height of the photon source is h = 0.1i?out above 
the disk plane. 

the high energy end is determined by the edge of the accretion disk. 

We have attempted to evaluate the outburst spectra in the 5-9 keV band again by 
implementing this Compton scattered He-like iron Ka line component in the model. In the 
fitting, the central energies of the iron lines are all fixed at their rest-frame energies (6.97 keV, 
6.67 keV, and 6.40 keV). The height of the plasma h is not constrained very well, and hence 
we also fix it at 0.1i?out; which is expected from the best-fit value of OP /2t: — 0.9. The result 
of the fit is shown in Fig. 13, and its best-fit parameters are summarized in table 6. Here 
we have employed a thermal bremsstrahlung as the continuum because it provides a better 
fit than a power law. The normalization of the Compton shoulder is tied to that of He-like 
iron Ka line according to the expectation from h = 0.1/2out- As a result of the fit, we found 
that the equivalent width of the broad 6.4 keV line is reduced from 2O8I27 ^'^ 1^0^33 eV. 
This value is, however, only marginally consistent with the equivalent width predicted from the 
given reflection scale (EW(fi/27r = 0.9i;^i) = 72i^^ eV). 

The obtained Gaussian a of the 6.4 keV line, on the other hand, is interpreted to 
represent the Keplerian velocity of the disk, which is 




(4) 
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Fig. 13. The outburst spectra in the 5-9 keV band with the best-fit model composed of a thermal 
bremsstrahlung and three iron lines including a Compton shoulder of the He-like Ka line in blue color. 

where Eq = 6.4 keV is the central energy of the fluorescent iron Ka line, i is the inchnation angle 
^37°. From the Gaussian a evaluated with the Compton shoulder model, cr = 0.073j;o;o57 keV 
(table 6), we obtain VK{r) = 5700^4400 km s~\ resulting in a radius of r = 3.512^/ x 10*^ cm for 
the I.IQMq white dwarf. The optically thin thermal plasma seems to exist on the disk within 
r ^ 7 X 10^ cm from the white dwarf. It may be intriguing to compare the spatial extension 
of the plasma thus obtained with that of an eclipsing dwarf nova in outburst or a nova-like 
variable. Based on the XMM-Newton observation of the nova-like variable UX UMa, Pratt et 
al. (2004) measured an elapsed phase for the the hard X-ray echpse transition to be A(p < 0.01. 
The binary parameters (Ritter & Kolb 2003 and references therein), on the other hand, indicate 
an orbital separation of a = 9.7 x 10^*^ cm. Hence, the plasma size is estimated to be < 10^ cm, 
in good agreement with the current estimation of SS Cyg in outburst. 

4-3. Elemental Abundances 

The elemental abundances of SS Cyg are summarized in table 3, which are plotted in 
the left panel of Fig. 14. All of the elemental abundances except for C are generally sub-solar. 
Done & Osborne (1997) determined the abundance of Si, S, and Fe separately, and that of the 
other elements in common by fltting ASCA outburst spectrum, which are also subsolar. We 
further resolved the abundances of Mg, Ne, O, N, and C. The right panel of Fig. 14, on the 
other hand, shows the abundances of 10 non-magnetic CVs obtained with the XMM-Newton 
observations (Pandel et al. 2005). These abundances are determined also by fits with the same 
model as ours. Their values are widely distributed from sub-solar to near solar abundances. 
The abundances of SS Cyg are within these distribution in general. 



22 



Table 6. Best-fit parameters of a thermal bremsstrahlung plus three iron Ka lines including a Compton shoulder to the 
outburst spectra in the 5-9 keV band. 



Component 


Parameter 


Best-fit value 


pon'l'inmTm 

\j\Jl± li 11J.LL LL1J.1 


rV-L I XvC V ( 


8 ()+"■'' 




n nvm f» 1 1 "7 f» 1 1 on T 


3 80+°-^*^ 

'^•°"-0.14 


H-like line 


enerffv fkeV) 


6.97 (fixed) 




normalization^ 


81+^-l? 

'-'•<JJ-_0.17 




EW (eV) 


78+11 

' "—16 


TTp-lilcp linp 




fi 67 ffixedl 




normalization^ 


4.93+°-^^ 

—0.30 




1—1 VV 1 C V 1 




Ocivn nton 


h/R f 


1 ffixedl 




EW (eV) 


42+2 


Neutral line 


enersrv fkeVl 


6.40 (fixed) 




a (keV) 


n7s+" 079 




normalization^ 






EW (eV) 


1401:^^ 




X'(d.o.f) 


205.7(157) 






1.31 



t In a unit of (3.02 x 10'^^ /{AirD)'^ EM), where EM is the emis- 
sion measure in a unit of cm~^. 
t In a unit of (xlQ-^ cm'^ keV"! s"!) at 1 keV. 
^ In units of (xlO~^ photons cm~^ 

5. Conclusion 

We have presented results of the Suzaku observations on the dwarf nova SS Cyg in 
quiescence and outburst in 2005 November. The X-ray spectra of SS Cyg are composed of 
a multi-temperature optically thin thermal plasma model with a maximum temperature of a 
few tens of keV, its reflection from the white dwarf surface and/or the accretion disk, and a 
6.4 keV neutral iron Ka line from the reflectors via fluorescence. High sensitivity of the HXD 
PIN detector and the high spectral resolution of the XIS enable us to disentangle degeneracy 
between the maximum temperature and the reflection parameters, and to determine the emis- 
sion parameters with unprecedented precision. The maximum temperature of the plasma in 
quiescence /cT^ax = 20.4+2^6 (stat.) ±3.0(sys.) keV is signiflcantly higher than that in outburst 
^^max ~ S.O+iJ keV. The elemental abundances of the plasma are close to the solar ones for 
the medium-Z elements (Si, S, Ar) whereas they declines both in lighter and heavier elements. 
Those of oxygen and iron are 0.46i[}:|]3 (stat.) ±0.01 (sys.)Z0 and 0.37l[]:[J3 (stat.) ±0.01 (sys.)Z0. 
The exception is carbon whose abundance is at least 2Zq even if we take into account all possible 
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Fig. 14. (left) Abundances of SS Cyg obtained with the Suzaku observations (§ 3.3). (right) Abundances 
of the 10 non-magnetic CVs obtained with the XMM-Newton observations in quiescence (Pandel et al. 
2005). Abundances are given relative to the solar values in Anders & Grevesse (1989) 

systematic errors. These trends are similar to other dwarf novae observed with XMM-Newton 
(Pandel et al. 2005). 

The solid angle of the reflector subtending over the optically thin thermal plasma is 
r2*°^/27r = 1.7 ± 0.2 (stat.) ±0.1 (sys.) in quiescence. Since even an infinite slab can subtend a 
solid angle of f2/27r = 1 over a radiation source above it, this large solid angle can be achieved 
only if the plasma views both the white dwarf and the accretion disk with substantial solid 
angles. Thanks to high energy resolution of the XIS, we have resolved a 6.4 keV iron Ka line 
into a narrow and broad components (significance of the broad component is ~99%), which 
also indicate contributions from both the white dwarf and the accretion disk to the reflected 
continuum spectra. The equivalent widths of them are both ~50 eV. Prom all these results, we 
consider the standard optically thin BL formed between the inner edge of the accretion disk and 
the white dwarf surface (Patterson & Raymond 1985) as the most plausible model to explain 
the observed large sohd angle. Prom the equivalent width of the narrow 6.4 keV component, 
the height of the BL from the white dwarf surface is h < 0.12i?wD- The total equivalent width 
of the 6.4 keV line (~100 eV) is consistent with that expected from Q^/27r, the iron abundance, 
and the incident illuminating continuum spectrum. 

The sohd angle of the reflector in outburst OP /2t^ — 0.9lo;4, on the other hand, is 
signiflcantly smaller than that in quiescence, and is consistent with an inflnite slab. Since the 
6.4 keV iron emission line is broad with no narrow component (^20% of the broad component), 
the reflection originates from the accretion disk. The accretion belt can also contribute to the 
reflection. The 6.4 keV line from the accretion belt is expected to be broad, which is consistent 
with the absence of the narrow 6.4 keV component. The EW of the 6.4 keV line is so large that 
it cannot be interpreted within a simple scheme of reflection from the disk. Even if Compton 
down-scattering of the observed He- like Ka line is taken into account, we can only flnd a 
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solution which marginally reconciles the large EW with the sohd angle of the reflector. We 
consider the optically thin thermal plasma in outburst as being distributed on the accretion 
disk. The Chandra HETG observation in outburst revealed that the He-like and H-hke emission 
hues from O, Ne, Mg, and Si are broad and their widths (~2000 km s~^) are consistent with 
those expected from the Keplerian velocity of the accretion disk (Okada et al. 2008). This 
fact suggests that the optically thin thermal plasma is anchored to the accretion disk and the 
accretion belt by magnetic field, for example, like solar coronae. 

Appendix 1. Low energy response of XISl 

In order to analyze spectra with XISl in the energy band ^0.5 keV where the energy 
response is affected by accumulation of contaminating material, we have checked the energy 
response of the XIS-1 with the PKS2 155-304 data taken between 2005 November 30 and 2005 
December 1 (seq.#700012010), which were taken close enough to our SS Cyg observations. 
Although the observation is originally planned to continue for a 60 ks effective exposure time, 
part of the observation suffered a light-leak accident and the CCD chip is irradiated by optical 
and UV photons^. After removing these time intervals, 39.4 ks data remain in total. We then 
created the same on-source and background regions as those adopted in the SS Cyg observations, 
and created the response file of XISl. The default correction was made for the contamination. 

Between 0.1 and 10 keV, PKS2155— 304 has been well represented by a curved spectrum 
with an energy slope gradually steepening from 1.1 to 1.6 (Giommi et al. 1998). We therefore 
fit our XISl spectrum with a broken power- law model 



where Fi and r2 are power law photon indexes, K is normalization factor, is a break energy, 
and see if the low energy calibration is enough for our analysis. The hydrogen column density 
to PKS2155 is 16.9 x 10^^ cm~^ from past infrared observations. The result of fit is shown in the 
left panel of Fig. 15. The fit below ~0.35 keV is not very well, and there is evidence of gain shift 
between 0.5-0.6 keV. Hence we apply an energy offset to the model by 1 eV step, and found 
that the — 11 eV shift gives the smallest reduced value. The remaining residuals are adjusted 
by introducing additional hydrogen column density and a carbon K edge at 0.2842 keV, which 
are both deemed associated with the contaminant of the XIS. The result of the fit is shown in 
the right panel of Fig. 15 and its best-fit parameters are summarized in table 7. The additional 
photoelectric absorption and the edge significantly improves the fit in 0.3-0.5 keV, and the 
spectrum down to 0.226 keV can be well represented by the broken power- law model. In the 
analysis of the Suzaku SS Cyg data, we always apply additional absorption column density of 
A^H = 8.2 X 10^^ cm"^ and the carbon edge with an optical depth of r = 0.88 at 0.2842 keV. 

^ http://wwwxray.ess.sci.osaka-u.ac.jp/ hayasida/ftp-files/XIS2/XIS_20060710b.pdf 



F{E) = K{E/l kcV)-^i {for E < E^) 

= KEl^-^' X {E/l keV)-^2 (/or E > E^), 



(Al) 
(A2) 
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Fig. 15. Spectral fitting of XIS-1 in the energy range of 0.1-10 keV before adjusting the energy scale 
(left) , and after adjusting the gain and applying additional absorption. After these treatments, the energy 
band above 0.226 keV can safely be used for spectral analysis. 

Table 7. Best-fit parameters of additional absorption component to fit the spectrum of PKS2155— 304. 



Model 


parameter 


Gain ofTsct (cV) 


-11 (fixed) 


T of carbon K edge 


0.88 ±0.05 


iVn (10^9 cm-2) 


8.2±0.7 


(d.o.f.) 


640 (547) 



The amount of gain offset is different from observation to observation. We have checked the 
offsets during the Suzaku observations of SS Cyg in quiescence and outburst in the same way 
as is described above to find that no offset 2 eV) is required for both observations. 
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